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Introduction : 





Déterminer la forme appropriée ou optimale des composants 
structuraux, qui soit à la fois de poids minimal et de rigidité maximale, est un 
problème de première importance pour l'ingénieur. Dans tous les domaines 
de la mécanique des structures, Гітрасі de la bonne conception d’une pièce 
est très important sur sa résistance, sa durée de vie et son utilisation en 
service. Des problèmes de ce genre se rencontrent fréquemment dans des 
domaines comme le génie civil, l'aéronautique, ou la construction 
automobile, pour lesquels toute économie de poids permet à la fois des 
économies de matière première mais surtout des gains très importants de 
performance (par exemple, une voiture légère consomme moins de 
carburant). En pratique, on demande aussi à la forme optimale de satisfaire 
d’autres contraintes technologiques (par exemple, aérodynamiques ou 
acoustiques) ou bien de faisabilité industrielle et de coût de fabrication. Le 
développement de l’art de l'ingénieur requiert des efforts considérables pour 
améliorer sans cesse les techniques de conception des structures. 
L’optimisation intervient de façon primordiale dans l'augmentation des 
performances et la réduction de masse des engins aérospatiaux et 
automobiles, entraînant ainsi de substantielles économies d'énergie. 


Habituellement, on procède par essais successifs, en testant des 
prototypes dont le design relève du savoir-faire et de l'intuition de 
l'ingénieur. Cette façon de faire manuelle est très coûteuse et imprécise. De 
plus en plus, elle est remplacée par des logiciels de modélisation numérique 
et d'optimisation, qui permettent d'analyser de nombreuses possibilités sans 
avoir à fabriquer de prototypes et qui automatisent la recherche de la forme 
optimale. Dans les codes de calcul usuels, cette automatisation est obtenue 
en représentant la forme par un nombre limité de paramètres descriptifs 
(généralement des points de contrôle sur les bords), et l'analyse des 
variations des performances par rapport à ces paramètres permet 
d'améliorer itérativement une forme initiale. 


Justification du choix : 





La topologie d’une structure étant une des premières décisions du processus de conception, on 
dispose, à ce moment, de peu d'informations ; on ne connaît souvent que certaines sollicitations 
prédominantes d’un côté et les fixations possibles de l’autre. Dès lors, dans de nombreuses 
applications, l'optimisation topologique a pour but de déterminer la nature et la connectivité des 
éléments constitutifs d’une structure dans un problème où seuls les conditions aux limites et le 
domaine spatial où il est possible de placer la structure, sont spécifiés. 





ідие : 


Problémat 





Cependant la méthode citée en introduction, dite d'analyse de sensibilité, 
présente deux inconvénients majeurs. Elle est en effet très coûteuse en 
temps de calcul car, si la forme s'éloigne trop de la géométrie initiale, il peut 
être nécessaire de modifier la discrétisation au cours des itérations. Par 
ailleurs, et c'est là leur plus grave défaut, le résultat obtenu dépend 
fortement du choix initial et de la finesse de la représentation. 


L'objectif qui nous intéresse est de trouver une forme optimale de la 
structure qui réalise le meilleur compromis entre sa résistance et son poids. 
Alors on se trouve face à la problématique suivante : 


Comment remédier aux problèmes de la méthode précédente en 
assurant une forme optimale d’une structure ? 





Solution proposée : 


On se propose donc de remédier aux inconvénients de la méthode d'analyse 


de sensibilité en proposant une nouvelle méthode dite : " Optimisation 
topologique par homogénéisation.”” Celle-ci permet d'optimiser 
automatiquement une structure sans restriction explicite ou implicite sur sa 


topologie, et évite les aléas numériques intrinsèques aux méthodes 
classiques, pour un coût de calcul très compétitif. 





L'optimisation est une branche des mathématiques cherchant à 
modéliser, à analyser et à résoudre analytiquement ou numériquement les 
problèmes qui consistent à minimiser ou maximiser une fonction sur un 
ensemble. Beaucoup de systèmes susceptibles d’être décrits par un modèle 
mathématique sont optimisés. La qualité des résultats et des prédictions 
dépend de la pertinence du modèle, du bon choix des variables que l'on 
cherche à optimiser, де l’efficacité de l’algorithme et des moyens pour 
le traitement numérique. 








Pour illustrer l’évolution des techniques d'optimisation structurale, on peut 
scinder arbitrairement l'optimisation des structures en trois grandes familles. 


Le dimensionnement automatique des 
structures : 

ne permet de modifier que la section droite 
ou l’épaisseur transversale des composants 
d’une structure dont la forme et la topologie 
sont fixées. Aucune modification du modèle 
géométrique n’est possible. 





l'optimisation de forme : admet des changements 

de forme compatibles avec une topologie fixée au préalable. 
L’optimisation de forme classique modifie 

la représentation paramétrique des frontières 

du domaine. En faisant bouger les frontières 

des domaines, on peut rechercher une meilleure 
solution parmi l’ensemble de toutes les structures 
obtenues par transformation homéomorphe de la 
structure originale. Dans ce cas, il est clair que l’on 
peut admettre un changement des dimensions 
transversales aussi bien qu’une modification de la configuration de la 
structure, mais il n’est certainement pas permis d’altérer la connectivité 
ou la nature des membres structuraux. 





l'optimisation topologique : permet de modifier plus fondamentalement 
la nature de la structure. Cette fois, la géométrie de 

la pièce est envisagée sans aucun a priori sur la 
connectivité des domaines ou des membres 
structuraux présents dans la solution. 
Optimiser la topologie conduit naturellement 
à déterminer d’une certaine manière la forme 
ou les dimensions transversales optimales de 
la structure, de sorte que certains auteurs lui 
attribuent aussi le nom d'optimisation de forme généralisée. 





L'optimisation de la topologie des structures constitue le sujet de ce rapport ! 








57 Définition : 

L'optimisation topologique est une méthode mathématique (et 
logicielle) qui permet de trouver la répartition de matière optimale dans 
un volume donné soumis à des contraintes. Elle se distingue notamment 
de l'optimisation de forme qui ne fait varier que la frontière de la pièce, 
qu'il faut avoir dessinée au départ. 

Les méthodes d'optimisation topologiques connues à ce jour 
permettent d'optimiser la résistance élastico-mécanique, la conductivité 
thermique ou certains problèmes d'écoulement fluide. 


7 Principe : 

Le principe de l’optimisation topologique est de considérer 
l’ensemble de tous les domaines possibles (inclus dans un domaine 
de référence Qo), sans restriction de régularité spatiale ou de 
connexité. Pour un tenseur l’élasticité L fixé, l’ensemble 2 est ainsi 
pris sous la forme 

2 = ЦІ, 0))0 с Qo}. 
Soit у la fonction caractéristique ди domaine Q, définie par 
x(x) = 1 si x E Q, et x(x) = 0 sinon. L'ensemble 2 peut être réécrit sous 
la forme : 


2 = {(L(x) = x(x)L, Qo) |X : Qo -> {0, 1}. 





У" Formulation 
La méthode dďd’homogénéisation minimise la fonction d'adaptation 
par rapport à de petites déformations sur le domaine de référence Q 
(dans lequel on cherche la forme optimale), en considérant les 
contraintes du problème. Cette méthode introduit des formes 
admissibles généralisées qui sont des matériaux composites 
caractérisés par deux paramètres : 
- la densité de matériau Ө(х) Є [0,1] 
- le tenseur homogénéisé qui représente la microstructure sous- 
jacente А"(х) 

Le problème d’optimisation relaxé ou homogénéisé (sans 
contrainte) s’écrit alors : min{ J (0, АТ) 
J Ла fonctionnelle à minimiser 

Cette méthode permet d'aboutir à un optimum global mais la 
solution obtenue peut être compliquée à exploiter pratiquement 
(matériau composite) 


Optimisation і topologlq ue 


1 


потодепеіѕг 
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Le problème le plus simple consiste à minimiser le poids d'une 
structure élastique, soumise à une ou plusieurs sollicitations, tout en 
respectant une contrainte mécanique (typiquement une borne 
maximale sur le travail des forces extérieures qui représente une 
mesure de rigidité globale). A l'intérieur d'un domaine de travail fixé, 
qui permet de tenir compte de contraintes d'encombrement, 
n'importe quelle forme est autorisée. Il n'y a donc aucune restriction 
sur les formes admissibles et le nombre de paramètres de design est 
illimité. En d'autres termes, cela revient à rechercher la répartition 
optimale de la matière et des trous dans ce domaine. 


Malheureusement, ce problème d'optimisation n'admet en général 
pas de solution si l'on n'élargit pas le cadre mathématique dans 
lequel il est naturellement formulé. D'un point de vue mécanique, 
cela s'explique par le fait qu'il est souvent plus avantageux de 
creuser, pour un même volume de matériau, un nombre infiniment 
grand de trous infiniment petits, plutôt que quelques gros trous. 
Dans les calculs numériques, ce phénomène se traduit sous la forme 
d'instabilités du résultat lorsque le maillage est modifié ou raffiné. 


Pour remédier à ces difficultés, il faut relaxer le problème en 
généralisant la notion de structure. Puisque les structures fines sont 
avantageuses, on s'autorise l'utilisation de matériaux composites 
obtenus par micro-perforation du matériau original, pour construire 
la forme optimale. 


Une forme ou structure généralisée est représentée, en chaque 
point, par une densité de matériau et par les propriétés mécaniques 
effectives du composite. Celles-ci dépendent fortement de la 
microstructure, c'est-à-dire de l'arrangement et de la forme des 
trous, ou pores, à l'échelle microscopique. La théorie de 
l'homogénéisation permet de calculer ces propriétés effectives et 





ensuite de déterminer la microstructure optimale, qui est, en chaque 
point, la plus rigide sous l'action d'un champ de contraintes donné. 


Cette relaxation par homogénéisation donne le bon cadre 
mathématique : on a transformé le problème de la recherche d'une 
structure composée de matériau et de trous, en celui de la recherche 
d'une structure composée en plus de matériau poreux. Les nouvelles 
inconnues sont, en chaque point, la densité de matériau poreux et sa 
microstructure. Dans la nouvelle formulation, la taille des trous n'est 
plus limitée par le pas du maillage. Bien sûr, la structure optimale 
peut comporter de larges zones composites de densité 
intermédiaire. Cependant, on sait comment la remplacer a posteriori 
par une forme classique (sans composite) aux performances proches 
et quasi-optimales pour obtenir ainsi un objet utilisable en pratique. 
Les exemples de l'utilisation de microstructures ne manquent pas 
dans la nature. Citons pour mémoire les os, dont la nature poreuse 
est adaptée aux efforts, ou les constructions en nid d'abeille, bien 
connues pour leur légèreté et leur solidité. 


L'utilisation de l'homogénéisation n'est pas seulement un artifice 
mathématique pour démontrer l'existence de formes optimales 
généralisées. Son importance va bien au-delà car elle conduit à une 
nouvelle classe d'algorithmes numériques, dits d'optimisation 
topologique, qui optimisent la répartition de matière sans aucune 
contrainte topologique. 









Modélisation 
Microstructure, Macroscopique 
є Réelle » Нотовепеіѕее 





Figure 1 : ргіпсірге de la méthode d’homogénéisation au sein d’un élément 


En figure1, la structure optimale réelle, souhaitée par l'algorithme 
d'optimisation, est modélisée par un élément homogène caractérisé 
par sa densité globale de matière. Le domaine d'optimisation va être 
composé d’un très grand nombre de ces éléments finis, ce qui 
représente autant de valeur de densité. L’algorithme d'optimisation 





a pour tâche de nous fournir Іа valeur де chacune de ces densités еп 
fonction des paramètres d'optimisation. 





L'algorithme que nous proposons se décompose en deux étapes : le 
calcul de la structure composite optimale, puis sa projection pour 
obtenir une forme classique. 


La première étape de calcul utilise un algorithme itératif : 
successivement, on résout un problème d'élasticité linéaire dans le 
domaine de travail, avec des paramètres de forme (densité, 
propriétés effectives des composites) donnés, puis, pour trouver les 
nouvelles valeurs de ces paramètres, on détermine la microstructure 
localement optimale à champ de contrainte donné, à l'aide des 
conditions d'optimalité. Comme l'étude de la formulation relaxée 
garantit l'amélioration de la performance à chaque itération, on 
répète cette séquence jusqu'à convergence. 


Cette méthode est stable et coûte peu cher en deux comme en trois 
dimensions d'espace. En effet, la mise à jour des paramètres de 
forme s'effectue localement grâce à des formules analytiques, et les 
calculs d'élasticité sont effectués par éléments finis, après 
discrétisation du domaine de travail. La forme et sa topologie sont 
donc capturées sur un maillage fixe, au contraire des méthodes de 
sensibilité pour lesquelles le maillage suit la forme à topologie fixée. 
La convergence s'obtient en quelques dizaines d'itérations et est 
indépendante des conditions initiales. 


La deuxième étape consiste à transformer la structure composite 
obtenue, en une forme classique, constituée uniquement de 
matériau pur. Pour ce faire, on pratique une pénalisation des 
composites : quelques itérations supplémentaires de l'algorithme de 
la première étape sont effectuées en forçant la densité à prendre des 
valeurs proches de 0 (trou) ou 1 (matériau pur). La forme finale ne 
contient plus de matériau composite, mais ses performances restent 
proches de l'optimum global donné par la première étape. Bien sûr, 
dans cette phase, la finesse du maillage conditionne la finesse des 
détails de la structure finale. 





En fait, nous concevons notre méthode comme un procédé 
fournissant une première ébauche de la forme optimale, qui définit 
sa topologie. On peut ensuite la coupler avec des méthodes 
classiques d'analyse de sensibilité pour raffiner localement les 
structures sélectionnées. Cette méthode est très versatile puisqu'elle 
permet d'optimiser une pièce soumise à un ou plusieurs 
chargements, et éventuellement avec des contraintes sur ses 
fréquences propres de vibration. 





Considérons un problème bidimensionnel inspiré de la conception 
d’un pont figure 2. On cherche à optimiser la forme d’un pont 
permettant de franchir une distance lo = 1.6 et supportant un 
chargement représenté par une force linéique f = - Є, appliquée еп 
у = 0.5, |x| < 1 (toutes les forces et longueurs sont exprimées sous 
forme adimensionnelle). Le volume occupé par le pont est fixé à 
0.25. On utilise la méthode d’homogénéisation pour chercher la 
forme de rigidité maximale vis-à-vis du chargement considéré. Pour 
ce faire, le domaine Qo est choisi égal au rectangle -1 < x < 1, 

0 < y < 1. Les conditions aux limites sont de type encastrement 
(déplacement nul) en 0.8 < |x| < 1, y = 0, et de type surface libre sur 
le reste de дО. 





Figure 2 : Recherche d’un pont optimal. Positionnement du problème 





у 0.3754 





Figure 3: Recherche d’un pont optimal. Solution relaxée р obtenue раг la méthode 
d'homogénéisation. 

La figure З montre la carte de densités p(x). Pour réaliser en pratique la 
distribution de matière représentée sur la figure 3, il faut imaginer une 
distribution de matériau composite (laminés), dont la microstructure évolue 
de façon continue en fonction de la position. Il est clair qu’une telle structure 
est difficile à réaliser en pratique. C’est pourquoi, en vue d’obtenir une forme 
facilement réalisable en pratique, un post-traitement est appliqué à la 
solution р(х) du problème relaxé, de façon à filtrer les densités 
intermédiaires et construire une distribution р" à valeurs dans {0, 1). 





0» 





Figure 4 : Recherche d’un pont optimal. Solution р" obtenue par homogénéisation, après 
filtrage 


Il est important de souligner que ce filtrage dégrade nécessairement la 
solution. De plus, la solution obtenue après filtrage est non unique et dépend 
notamment du maillage utilisé. 
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Perspectives d’élargissement 


Le développement constant des techniques de conception assistée par 
ordinateur et des stratégies d'optimisation s'inscrit dans ce cadre. 
L’optimisation des structures soulève depuis plus de vingt ans le plus мії 
intérêt. Encore trop peu appliquée aux techniques classiques de bureau 
d’études, elle s’y intègre progressivement au fur et à mesure que s'accroît sa 
fiabilité. Parti des problèmes les plus simples, le champ d’application de 
l'optimisation structurale s'étend aujourd’hui à de nouveaux défis toujours 
plus intéressants. L’optimisation de la topologie des structures est 
maintenant un sujet en pleine ébullition. L'établissement de ses bases 
théoriques remonte au milieu de la dernière décennie, tandis que les 
premières applications pratiques ne datent que de la fin des années 80. 
Parallèlement au développement de nouvelles microstructures et à 
l'émergence d’approches alternatives, le champ d’application de 
l'optimisation topologique, fondé sur la méthode d’homogénéisation, n’a pas 
cessé de s'étendre. Initialement limitée aux membranes et à une seule 
sollicitation statique, la méthode d'homogénéisation a été généralisée а 
d’autres types de structure ainsi qu’à d’autres types d’analyse et à d’autres 
conditions de conception. L'outil d'optimisation topologique par 
homogénéisation est maintenant arrivé à une maturité suffisante pour que 
l’on puisse voir les premières applications industrielles. Le problème de la 
poutre dite "МВВ" est certainement une des premières applications de taille 
réaliste. A l'heure actuelle, ce sont principalement le domaine aéronautique 
et plus encore le secteur de l'automobile qui semblent être les principaux 
champs d’application de l'optimisation topologique. 


Application au Maroc : 





Au Maroc, comme dans le monde entier, l’optimisation topologique par 
homogénéisation est une technique qui а connu une grande évolution 
pendant les dernières années puisqu’elle met en œuvre des méthodes 
efficaces qui ont pour but l'optimisation des structures afin de minimiser les 
coûts de production et l'énergie. 

Cette méthode est très utilisée dans les différents secteurs d'activités et plus 
fréquemment en mécanique (matériaux de construction), aéronautique et 
l'automobile, qui sont parmi les piliers de l’économie marocaine. 





CONCLUSION : 





L'approche topologique est le plus grand défi soumis à l'optimisation de 
structure. La définition topologique d’une pièce affecte en premier lieu ses 
performances, elle est donc du premier ordre en terme d’optimalité. En effet 
les résultats ne seront jamais optimaux si la topologie initiale n’est pas la 
bonne. Les méthodes d'optimisation topologiques actuelles ne sont pas 
capables de définir définitivement le forme d’une pièce, mais elles 
représentent de très fiables préprocesseurs pour les autres méthodes plus 
précises et plus correctes du point de vue conception technologiques. 

L'optimisation topologique est souvent perçue comme un outil de 
conception automatique qui va permettre principalement d’alléger ипе 
pièce. Cependant, sans l'expertise et le savoir-faire du concepteur, Іа 
conception topologique peut aboutir à des conceptions optimisées 
mécaniquement mais difficiles à réaliser. L’optimisation topologique, si elle 
peut être utilisée pour des réductions de masse, peut aussi être utilisée à bon 
escient dans le cadre d’une reconception classique et aboutir à des 
géométries optimisées inédites, en rupture avec la conception initiale. 
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